Langevin dynamics in crossed magnetic and electric fields: Hall and diamagnetic 

fluctuations 
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Based on the classical Langevin equation, we have re-visited the problem of orbital motion of a 
charged particle in two dimensions for a normal magnetic field crossed with or without an in-plane 
electric bias. We are led to two interesting fluctuation effects: First, we obtain not only a longitudinal 
"work-fluctuation" relation as expected for a barotropic type system, but also a transverse work- 
fluctuation relation perpendicular to the electric bias. This "Hall fluctuation" involves the product of 
the electric and the magnetic fields. And second, for the case of harmonic confinement without bias, 
the calculated probability density for the orbital magnetic moment gives non-zero even moments, 
not derivable as field derivatives of the classical free energy. 
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Consider the classical motion of a charged particle in a 
plane under the influence of a normal (out-of-plane) mag- 
netic field crossed in general with a parallel "in-plane" 
electric field, in the presence of a dissipative coupling to 
the environment (bath) . This 2D motion has two notable 
general features - first, the cyclotron circular orbits about 
the magnetic field under the Lorentz force that classically 
does no work, and second, there is the well known trans- 
verse drift identified as the Hall effect which is a dis- 
sipative transport (non-equilibrium) phenomenon. The 
orbital motion holds surprises— the best known of these 
being the absence of classical orbital diamagnetism in 
equilibrium, as embodied in the classic theorem due to 
Bohr and van Leeuwen, and, indeed, regarded by some 
to be a surprise of theoretical physics A notable 
point here is the subtle role played by the spatial con- 
finement (the boundary) that interrupts the otherwise 
complete circular (diamagnetic) orbits leading to a skip- 
ping cuspidal motion at the boundary which is retrograde 
(paramagnetic) and exactly cancels the bulk diamagnetic 
contribution 0, H ,_!)]• Diamagnetism, or rather its ab- 
sence classically, |5| is, of course, a well known equilib- 
rium phenomenon. We, however, find interesting fluctu- 
ation effects here. Thus, while the orbital diamagnetic 
moment is zero in the mean, its higher (even) moments 
do exist under confinement, but can not be obtained as 
field derivatives of the classical free energy. The other 
surprising effect manifests in the non-equilibrium steady 
state as a transverse work-fluctuation relation in addi- 
tion to the well known longitudinal work-fluctuation the- 
oretic results of much current interest 
the light of these, we re-visit the classical problem of 
orbital motion of a charged particle in two dimensions 
in the presence of a normal magnetic field. First, we 
will consider the work-fluctuations in the presence of a 
parallel (in-plane) electric bias. Keeping in mind the 
non-equilibrium steady-state situation of interest here, 
we follow the Einsteinian statistical approach based on 
the Langevin equation whose solution in the long-time 



limit gives the non-equilibrium steady-state results for 
the system driven externally (e.g., biased electrically). 
The main results are the following: For non-zero elec- 
tric bias, a "work fluctuation" expression obtains as, of 
course, expected for a barotropic-type system. In addi- 
tion, however, a transverse work-fluctuation is also ob- 
tained perpendicular to the electric bias - a "Hall fluc- 
tuation" theorem. This "Hall fluctuation" involves the 
product of the electric and the magnetic fields. Inasmuch 
as classically the magnetic field does no work on a mov- 
ing charge, this cross effect (Hall fluctuation) seems to be 
related to the time-reversal symmetry breaking effect of 
the externally applied magnetic field. This is quite dis- 
tinct from the time-reversal asymmetry (irreversibility) 
arising out of dissipation. It is apt to point out at this 
stage that the physical scheme we have in mind here is 
experimentally an obvious variation on the well known 
Haynes-Shockley [l(| set-up for studying the drift con- 
comitant with the spread (diffusion) of minority charge 
carriers photo-injected in an electrically biased semicon- 
ducting sample. Indeed, the fact that the minority car- 
riers can disappear through the electron-hole recombina- 
tion process, provides for a study of finite life-time effects 
in our non-equilibrium steady state This should be 
an interesting way of probing the effect of confinement, or 
the boundary, that may not be effective when the carrier 
life-time is very short. 



Fluctuation theorem - The classical Hall bar : Con- 
sider a classical 2D system of noninteracting electrons 
(charge — e) under the externally applied crossed electric 
and magnetic fields, in the presence of a dissipative envi- 
ronment at temperature T (the bath) in the high friction 
limit. Let the electric field (E) be in the x— direction 
and the magnetic field (B) in the z— direction. We can 
now write down the Langevin equations of motion for 
the electron in the overdamped limit, i.e., ignoring the 
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inertial effects 11[ as: 



73; = -eE y + Vx(t) 

c 

, • eB ■ 

and 7?/ = —-x + rj y {t) 

where 7 is the friction coefficient due to the underlying 
bath and r\ x (t) and rj y (t) are the concomitant random 
(stochastic) forces generated by the bath. The electrical 
drift mobility is defined as e/7 for electrons. We assume 
as usual the random forces to be Gaussian white noise 
with (ri a ,p(t)) = and (v a (t)^)) = vl S(t - f ) S af3 . 
Here (...) denotes average over the stochastic ensemble at 
time t. Solving the above two linearly coupled equations 
for the velocities x and y, we obtain 
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with eB/jc = and eE/j = a. In the overdamped 
limit, the information of the electron's motion is cap- 
tured in the probability distribution P(x,y,t) of the x- 
and y-positions at time t. We are interested in finding 
the marginal probability distributions (P(x,t), P(y,t)) 
of the displacements x and y respectively: P(x, t) = 
J P(x,y,t)dy and P(y, t) = J P(x,y,t)dx. From the 
van Kampen lemma [12j, we can write down a conti- 
nuity equation for the marginal density ir(x,t) evolving 
stochastically under Eq.{T]) with P(x,t) = (Tr(x,t)) as 



dir(x, t) 
dt 
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The Fokker-Planck equation for the marginal probability 
distribution P{x, t) is now obtained by the noise aver- 
aging of Eq.([3]). We apply Novikov theorem [l2| for the 
Gaussian noise and find 



dP(x,t) a dP(x,t) 
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dt 
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(4) 
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The last equation describes the diffusion of electrons 
with drift. Similarly, the Fokker-Planck equation for the 
marginal probability distribution P(y,t): 



dP(y,t) 
dt 



aP dP(y,t) 
1 + p 2 dx 
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d 2 P(y,t) 



2 7 2 (1 + p 2 ) dy 2 



It is easier to solve the differential Eqs. (|4|5p by the 
Fourier transform method with initial conditions P(x, t = 
0) = 5(x) and P(y,t = 0) = S(y). It is given as 



P(x,t) 
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The marginal probability density P(x, t) is a Gaussian 
distribution with mean (x) = —at/ (1 + B 2 ) and variance 
((x 2 ) - (x) 2 ) = v 2 th 2 (l + B 2 ). Here v xd = -a/(l + B 2 ) 
is the drift velocity along the a;— direction. It can also be 
derived from Eq. |T]) by taking ensemble average on both 
sides. The above solution, of course, reduces to ordinary 
diffusion in the limit E = and B = with the identifi- 
cation r] 2 = 2D7 2 , where D is the diffusion constant. We 
can now readily derive the usual barotropic-type work 
fluctuation relation in the direction of electric field: 



P(-x,t) 
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with the Einstein relation, Dj = ksT. The interpreta- 
tion of the above relation is straightforward. Though the 
system is always far from thermodynamic equilibrium (as 
the electron continuously dissipates energy in the envi- 
ronment), it reaches a mechanical equilibrium asymptot- 
ically under the combined effect (forcing) of the electro- 
magnetic fields and the viscous drag acting in opposition. 
After the system attains the steady state, if the electron 
is at a position x\ , say, at t = 0, then P(x, t) is the condi- 
tional probability of finding the particle at position X2 at 
later time t with the displacement x = xi — x\. Now, we 
reverse the positions keeping time direction as before and 
start at t = from the position X2, then P(—x, t) is the 
conditional probability of finding the particle at position 
x± at time t. The main feature of the relation ^ is that 
right-hand side is time independent. The relation ^ can 
also be derived from the eqilibrium probability distribu- 
tion (depicting microscopic reversibility) [l3| . Using now 
the normalisation condition of the probability distribu- 
tion, we obtain another useful relation from Eq.®: 
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P(-x,t)dx = l (7) 
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The last two relations (Eq.((6|),Eq.([7])) are similar to the 
generalised fluctuation-dissipation theorems which were 
derived in a very different context in Ref . [14| . 

Next, we investigate the rather novel and interesting 
transverse fluctuations. Again, using the Fourier trans- 
form method, we solve Eq.([5]) for the marginal probabil- 
ity distribution (P(y,t)) of y giving 
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where the drift velocity along the y— direction is v yc i — 
—a(3/(l + B 2 ). The ratio of the probability densities for 
the transverse fluctuations is given by 
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FIG. 1: (color online). Plot of the probability density P(M) 
of Eg . (llOfl vs. magnetic moment M for different values of if. 
Here M is scaled by the Bohr magneton (/is). 



where the interpretation of P(y,t) and P(—y,t) are sim- 
ilar to that for the x-direction. We call the last relation 
"Hall fluctuation" theorem -it involves the product of 
the crossed electric and magnetic fields. The magnetic 
field, of course, does no work on a charge moving in the 
xy-plane. This cross "Hall fluctuation" effect can be re- 
lated ultimately to the time-reversal symmetry breaking 
effect of the applied magnetic field, which is different from 
the time-irreversibility introduced through the dissipa- 
tion (7). We can interpret it variously as arriving from 
the transverse Hall voltage (the Lorentz force), or effec- 
tively as a magnetoresistance in the Hall geometry which 
occurs due to the enhancement of the path length of elec- 
tron's motion caused by the magnetic field. Finally, we 
give a generalised fluctuation-dissipation theorem-like ex- 
pression in the presence of the magnetic field: 



Diamagnetic fluctuations : Next, we turn to the dia- 
magnetic fluctuations. We derive and briefly discuss the 
classical equilibrium fluctuations of the orbital diamag- 
netic moment (which is known to vanish identically in 
the mean). More specifically, we find the equilibrium 
probability density for the orbital magnetic moment of 
a charged particle in two dimensions for a normal mag- 
netic field and a harmonic confinement. The probability 
density, when properly scaled, turns out to be universal 
and peaked about a zero mean value. Here, however, we 
must retain the inertial effect (non-zero electron mass). 

Consider the orbital motion of the electron in the xy- 
plane for the normal external uniform magnetic field B 
along z— directions and a harmonic confinement in the 



xy-plane of strength fco- The Hamiltonian 
1 , eBy^ 1 
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where we have used the symmetric Landau gauge (A = 
B x r/2) for the vector potential A. The orbital diamag- 
netic momemt can be expressed as 



M(t) 
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We are interested in the asymptotic (t — > 00) distribution 
of M(t). Now, we can evaluate the probability density 
P(M) of classical M in the limit t — > 00 (i.e. in equilib- 
rium) through the usual method of finding the equilib- 
rium distribution at finite temperature T: 
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Here Z is a normalisation constant to be determined 
through P(M)dM = 1 (infact Z is the equilibrium 
partition function for this model). After some simple al- 
gebra we obtain the probability density 



P(M) 
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where fi B = eh/2mc is the Bohr magneton and ujq — 
y/ko/m. As expected (M) — 0, i.e., orbital diamag- 
netism is identically zero in a confined system. But, 
we do find finite thermal fluctuations of the orbital dia- 
magnetic moment in confinement. In FIG. [1] we plot 
P{M) for different values of the dimentionless param- 
eter ip — fuoo/k B T. It is interesting to note that the 
classical partition function (Z), and hence the free en- 
ergy (— k B T\nZ) are independent of the magnetic field 
B giving vanishing orbital diamagnetism in the mean. 
But the moments of the orbital magnetic fluctuation of 
even order are all non-zero. Clearly, then these orbital 
magnetic moments can not be derived through the usual 
field derivatives of the classical free energy as would be 
the case for "permanent magnetic moments" intrinsic to 
the particles. The classical equilibrium simulation of the 
above model with the Langevin heat bath can be shown 
to be consistent with the above thermal fluctuations of 
the magnetic moments. We simulate the coupled set of 
Langevin equations 

. eB . , 

mx — —"fx y — kqx + rixit) 

c 

. eB . , 

and my = —7?) H x — hoy + rj y {t) 

where again 7 and 77 are respectively the friction and the 
noise, related through {ila(t)r]0(t')) = 2jk B TS(t—t r ) S a p. 
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FIG. 2: (color online). Plot of the probability density P(M) 
vs. magnetic moment M obtained in simulation for same 
values of p as in FIG. [1] Peak of the curves decreases with 
decreasing (p. We have used y = 1, m = 1 and again M is 
scaled by fiB- Inset P(M) for ip — 1.0; and for different values 
of 7 and B showing no effect. 



Here we use the velocity- Verlet algorithm for the time 
evolution of the above equations and find steady-state 
equilibrium distribution of the magnetic moment given 
by Eq.©. We plot the distribution P(M) in FIG. for 
the same parameter values of <p as in FIG. [1] Also, we 
confirm through our simulation that the P(M) is inde- 
pendent of 7 and B (see inset of FIG. [2]). 

Concluding remarks : In this paper we have re-derived 
a barotropic-type work-fluctuation relation along with a 
new transverse fluctuation relation for the case of the 
classical motion of a charged particle in static homoge- 
neous crossed magnetic and electric fields in the presence 
of dissipation. This is interesting inasmuch as classically 
the magnetic field does no work on a moving charge par- 
ticle. Recently, there have been some studies of fluctua- 
tion theorems [15| | in time- varying electromagnetic fields. 
But, our approach and the context are quite different 
from these. As our treatment is based on the classi- 
cal Langevin equations involving stochastic fluctuating 
forces and the concomitant dissipation that neglect quan- 
tum statistics and the band structure effects [lfj, it is 



expected to be appropriate for a material system which 
is electronically non-degenerate, i.e., has a low carrier 
density at a relatively high temperature, and has a low 
carrier mobility. Clearly, a polar semiconductor with 
strong electron-phonon interaction is indicated. The in- 
ertial mass occurring in our Langevin equations is, of 
course, to be replaced by the effective mass relevant to 
the bottom (top) of the conduction (valence) band for the 
electron (hole) as the band is expected to be close to be- 
ing parabolic there. We have also derived the probability 
density for the classical diamagnetic moment giving non- 
zero orbital magnetic moments of even order. The latter 
can not be obtained as field derivatives of the equilibrium 
free energy which is clasically field independent. 

One of us (D.R.) thanks A. Dhar and R. Marathe for 
fruitful discussions. 
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